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Abstract 

It is shown that, in the scahng regime, transport properties of quantum wires 
with off-diagonal disorder are described by a family of scaling equations that 
depend on two parameters: the mean free path and an additional continuous 
parameter. The existing scaling equation for quantum wires with off-diagonal 
disorder [Brouwer et al., Phys. Rev. Lett. 81, 862 (1998)] is a special point in 
this family. Both parameters depend on the details of the microscopic model. 
Since there are two parameters involved, instead of only one, localization in a 
wire with off-diagonal disorder is not universal. We take a geometric point of 
view and show that this nonuniversality follows from the fact that the group of 
transfer matrices is not semi-simple. Our results are illustrated with numeri- 
cal simulations for a tight-binding model with random hopping amplitudes. 

PACS numbers: 72.15.Rn, 72.10.Bg, 11.30.Rd 
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Universality is a key concept in any approach to study localization in disordered systems. 
In the scaling theory of localization, it is commonly believed that the statistical distributions 
of the conductance, or of energy levels and wave functions, are entirely determined by the 
fundamental symmetries and the dimensionality of the sample [|l|. Once symmetry and 
dimensionality are taken into account, all microscopic details of a sample can be represented 
by a single length scale the "mean free path", such that on length scales L larger than 
the sample is completely characterized by the ratio L/^. The concept of universality is 
the cornerstone of various field-theoretic, diagrammatic, and random-matrix approaches to 
localization 

For a quasi-one dimensional geometry, i.e., for quantum wires, and for weak disorder 
(mean free path I is much larger than the Fermi wavelength A) the statistical distribution 
of the conductance can be described by the transfer matrix approach of Dorokhov , and 
Mello, Pereyra, and Kumar (DMPK). In this approach, the transport properties of the 
quantum wire are described in terms of a scaling equation for its transmission eigenvalues, 
the so-called DMPK equation 0. Hiiffmann and Caselle ^ have provided the DMPK 
equation with a geometric foundation by reformulating it in terms of a Brownian motion of 
the transfer matrix on a symmetric space, a certain curved manifold from the theory of Lie 
groups and Lie algebras The existence of a unique natural mathematical framework to 
describe Brownian motion on symmetric spaces provides the geometric counterpart of the 
observed universality of the localization properties of disordered quantum wires. 

In a recent work, two of the authors, together with Simons and Altland [0, have pro- 
posed an extension of the scaling approach of DMPK to quantum wires with off-diagonal 
disorder (e.g. a lattice model with random hopping amplitudes and no on-site disorder). At 
the band center e = such a system has an extra chiral or sublattice symmetry that is not 
present in the standard case of a wire with diagonal (potential) disorder. Therefore they 
belong to a different symmetry class, which is referred to as the chiral symmetry class. Chi- 
ral symmetry also plays an important role for two-dimensional Dirac fermions in a random 



vector potential [0-h3|, the lattice random flux model |14-16|, non-Hermitian quantum 



mechanics [|T^|T8[, supersymmetric quantum mechanics |T^, diffusion in a random medium 
pO| , and in certain problems in QCD ||21|| . 

In Ref. |T^, the extension of the DMPK equation to the chiral symmetry class was de- 
rived from a simple microscopic model. Here we discuss its geometric origin. This proves 
to be an exercise with implications that reach far beyond the construction of a mere chiral 
parallel to the geometric foundation of the DMPK equation of Refs. Upon inspec- 

tion of the geometric structure underlying the "chiral DMPK equation", we find that the 
localization properties of a disordered wire with off-diagonal disorder are not universal; the 
geometric approach allows for a one-parameter family of scaling equations for the transmis- 
sion eigenvalues. The scaling equation that was originally found in Ref. [0 is a special point 
in this family. Below we discuss the geometric approach in more detail, identify the origin 
of this non-universality, and illustrate the results with numerical simulations of different 
microscopic models. 

We start with a brief summary of the ideas that lead to the standard DMPK equation and 
its geometric interpretation. The cornerstone of this approach are the symmetry properties 
of the transfer matrix M of a disordered wire. For definiteness, we focus on the lattice 
Anderson model in a geometry of coupled chains, see Fig. |I|. In this model, the 
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FIG. 1. Lattice model for a disordered sample consisting of = 3 coupled chains (dot- 
ted region) attached to two ideal leads. In the two leads, the wave function is represented by 
two A^-component vectors of amplitudes for the left and right moving waves respectively. In the 
disordered region, the Hamiltonian contains a random on-site potential. 



Hamiltonian consists of nearest neighbor hopping terms and a random on-site potential. 
We restrict our attention to spinless particles. The wire consists of a disordered region and 
of two ideal leads consisting of N (uncoupled) chains. In the leads, the wave function is 
represented by an A^-component vector ip± for the amplitudes of left (+) and right (— ) 
moving waves. The wave functions to the left and right sides of the disordered sample are 
related by the 2 x 2N transfer matrix M , 

C-) right ^(v'-)ieft' 

In this basis of left and right movers, flux conservation implies that the transfer matrix M 
obeys 

M^EaM = S3, (2) 

where S3 = as ® Hat, (73 being the Pauli matrix and Hat the N x N unit matrix. We 
distinguish between the presence and absence of time-reversal symmetry, labeled by the 
symmetry parameter (3 = 1,2, respectively. In the presence of time-reversal symmetry, M 
further satisfies 



SiM*Si = M, (3) 
where Si = cti ® Hat. The transfer matrix M can be parametrized as |^,22 



_ I U \ / cosh a; sinhx \ f ^ ^ \ /a\ 



U' J \ sinh X cosh x J \ V' 

where U, U', V, and V are unitary matrices, and x is a diagonal matrix with diagonal 
elements Xj. For (3 = 1, V = V* and U' = U*. The unitary matrices U, U', V, and V 
serve as "angular coordinates" for M, the parameters Xj serve as "radial coordinates" [the 
eigenvalues of MM^ are exp(±2a;j)]. The radial coordinates Xj are related to the transmis- 
sion eigenvalues Tj by Tj = 1/ cosh^Xj. They determine the dimensionless conductance g 
through the Landauer formula, 
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N N ^ 

^ ^ cosh Xj 

In the original derivations of the DMPK equation the wire is divided into thin 

shces, and the transfer matrix M of the entire wire is found by multiphcation of the transfer 
matrices of the individual slices. For weak disorder, the parameters Xj of the transfer matrix 
undergo only small changes xj — > xj + 6xj upon each such multiplication. One can view 
this process as a "Brownian motion" for the parameters xj, the length L of the wire serving 
as a fictitious time. With the choice of a maximum information entropy distribution for 
the transfer matrix of the slice one can write down the corresponding Fokker-Planck 
equation, which has the form 

dL ^ dxj dxj 
J = Y\ sinh(2xj) Y\ sinh'^(a;j — Xk) sinh^(xj + Xk). (6) 

The proportionality constant D is determined by the details of the microscopic model.Q 

Hiiffmann pointed out that there is a beautiful geometric structure underlying the scaling 
equation (^. He observed that Eqs. (^-(^ express that the transfer matrix M is a member 
of a Lie group Gp, with[] 

Gi ~ X Sp{N, M), ^2 ~ U{N, N) ^ f/(l) X SU{N, N), (7) 

and proposed to describe the L evolution of M as a "Brownian motion" on the manifold Gjj 
0. For a rigorous formulation of this Brownian motion process, two further steps have to 
be taken. We would like to repeat them here, as we need to reconsider them when we deal 
with the case of a quantum wire with off-diagonal disorder. 

1. The Lie groups are not semi-simple: they are the direct product of the two com- 
ponents Z2 and Sp {N, M), or f/(l) and SU{N, iV), for = 1 or 2 respectively. For the 
product MM^ , and hence for the radial coordinates xj, only the non-compact compo- 
nents SplNjWj and SU{N,N) are relevant, so that we may restrict our attention to 
the semi-simple Lie groups Sp{N,M.) and SU{N, N). The remaining components Z2 
and f/(l) correspond to the sign or phase of detM and do not affect the xj. 

2. The algebraic structure on the semi-simple Lie groups Sp{N,M.) and SU{N,N) gives 
rise to a natural metric 0. However, since this natural metric is not positive def- 
inite, it cannot be used to define a Brownian motion process. The problem is 



^ In the DMPK equation, D^^ = 2{/3N + 2- (3)£, where £ is the mean free path, see e.g. Ref. 

^ U{N,N) is the group of complex matrices M with M^Ti^M = S3; SU{N,N) is its subgroup 
of matrices with unit determinant; Sp{N, M) is the group of real 2A'" x 2N matrices M that obey 
M'^EsM = S2, and Z2 = {±1}. To see that Gi ~ Z2 x Sp (iV,]R), see e.g. Ref. @. 
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solved by dividing out the maximal compact subgroups U{N) or S{U{N) x U{N)) 
for /3 = 1 or 2, respectively. The resulting coset spaces Si = Sp{N,'R)/U{N) and 
5*2 = SU{N, N)/ S{U{N) X U{N)) are called symmetric spaces and have a natural 
positive definite metric 0. In the parametrization (H), the procedure of dividing out 
the subgroups U{N) or S{U (N) x U {N)) corresponds to the identification of all transfer 
matrices M, M' for which the product M~^M' is of the form 

M-'M' = ( °, ) , (8) 

where V and V are arbitrary unitary matrices {V = V* for (3 = 1 and det V'V = 1 for 
(3 = 2). In shifting to the symmetric spaces Sf3 no information on the radial coordinates 
Xj is lost, because they are well-defined in each equivalence class. 

The symmetric spaces 5*^3 admit a spherical coordinate system whose radial coordinates 
Xi, . . . ,X]sf are equal to the radial coordinates xi, . . . ,xn of the transfer matrix M, 
cf. Eq. (^. The Brownian motion of the xj is described by the radial part of the 
Laplace-Beltrami operator on Sp. This radial part is known from the literature f^-Q, 
and is found to be identical to the differential operator in Eq. (||). As a result, we 
find that their probability distribution P{xi, . . . ,xj\f; L) satisfies the Fokker-Planck 
equation (P), where the constant D is the diffusion constant on Sp. The appearance 
of a single diffusion constant D signifies the universality of the localization properties 
in disordered quantum wires. 

Let us now consider a quantum wire with off-diagonal disorder (random hopping), at the 
band center e = 0. We consider the same geometry as in Fig. the only difference being 
that now the randomness is in the hopping amplitudes between neighboring lattice site, the 
on-site potential being zero everywhere. On a bipartite lattice (which is the case we consider 
here, cf. Fig. 0), and with only off-diagonal (hopping) randomness, localization is different 
from the standard case described above because of the existence of an additional symmetry, 



known sublattice or chiral symmetry ||10| 



SiMSi = M. (9) 

Hence the unitary matrices in the parametrization (^ satisfy U = U' and V = V . Let us 
now consider the extension of the DMPK equation (^) that includes the chiral symmetry 
(y), following Hiiffmann's geometric approach. 

The Lie groups of transfer matrices M that obey the chiral symmetry (^ areQ 

Gf ~ GL{N, M) ~ ^2 X M+ X SL{N, M), 
Gf - GL{N, C) - U{1) X M+ X SL{N, C). 



^ GL[N, M) and GL{N, C) are the multiplicative groups N x N matrices with real and complex 
elements, respectively; SL{N,E.) and SL{N,C) are their subgroups of matrices with unit deter- 
minant; is the multiplicative group of the positive real numbers. The origin of these transfer 
matrix groups is explained below Eq. ([17|). 
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The construction of a Brownian motion process for M on proceeds with the same two 
steps as in the standard case of diagonal disorder. However, there is an important difference. 
Unhke the transfer matrix groups Gp discussed above, the transfer matrix groups for off- 
diagonal disorder contain two non-compact factors, and SL{N, M), or and SL[N, C), 
for /5 = 1 or 2, respectively. (This difference was noted by Zirnbauer in a field-theoretic 
context ||2^.) Both of them determine the radial coordinates Xj: The factor M"*" describes 
the position of the average x = [xi -|- . . . -|- xn)/N, while the special linear group SL{N) 
is connected to the relative positions of the Xj. One therefore has to consider two different 
Brownian motion processes: One on M+, to describe the L evolution of x, and one on the 
symmetric space Sf" = SL{N,R)/SO{N) or = SL{N,C)/SU{N), to describe the L 
evolution of the differences of the radial coordinates xj. [The symmetric spaces S"^ are 
obtained after dividing out the maximal compact subgroup of SL{N, M) or SL{N, C), as in 
the standard case.] Since these two Brownian motion processes have two different and a priori 
unrelated diffusion constants, one needs their ratio as an extra parameter to characterize 
the distribution of the transfer matrix. This is the absence of universality in quantum wires 
with off-diagonal disorder that we announced in the introduction. 

In this most general case, there is a one-parameter family of Fokker-Planck equations for 
the distribution P{xi, . . . ,xn', L) of the radial coordinates Xj. This family of Fokker-Planck 
equations is constructed in two steps. First, Brownian motion on M"*" results in a simple 
diffusion equation for the distribution Px{x; L) of the average x, 

dPs, _Dr d"^ 

Here Dj^/N is the diffusion coefficient for the Brownian motion on M"*". (We chose to write 
the diffusion coefficient as Dji/N for later convenience; it is the diffusion coefficient that 
one would obtain for the average x if the N variables Xj would diffuse independently with 
diffusion coefficient Dr.) Second, Brownian motion on S'^^ is described in terms of radial 
coordinates y^, j = 1, . . . , A^, with the constraint ?/i + . . . + ?/Ar = [1,0. They correspond to 
the radial coordinates Xj of the transfer matrix M Using the explicit form 

for the Laplace-Beltrami operator on the symmetric spaces S'^', the Fokker-Planck equation 
for the probability distribution Py{yi, . . . ,yN; L) of the yj reads 

f)P ^ f) f) 

j=l j<k 

where we have to restrict to the subspace yi + . . . + y^ = 0. Hence, for the distribution 
P{xi, . . . ,xn', L) of the radial coordinates Xj = x + yj of the transfer matrix M we find the 
Fokker-Planck equations 

W-^'^^ar, [''^ ' ^ n - -A (12) 

j,k — 1 j'^k 

where D^h = Ds and r] = Dr/Ds- Equation (|I2|) contains the full one-parameter family 
of scaling equations that describes locahzation at the band center e = in quantum wires 
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with off-diagonal disorder. The case 77 = 1 (i.e., equal diffusion constants Dr and Ds) 



corresponds to the equation that was derived before in Ref. |]10 



A solution of the scaling equation (|T2]) in the localized regime L ^ can be obtained 
by standard methods 0. The distribution of the radial coordinates is Gaussian, with average 
and variance given by 



{xj) = {N + 1- 2j)/3LDch, {xkXj) - {xj){xk) = 2 Sjk - 



V 



N 

For even A^, the distribution of the conductance g is log-normal, with 

1/2 



ch- 



(13) 



{lng) = -(3s + 2 



i 1-2!^ 

TT V N 



var Ing = 2 



1+1- 



TT 



1 -2 



1 — rj 
N 



s + 0{l), 



(14) 



where s = 2LDch- For odd A^, the fluctuations of InG are of the same order as the average, 

1/2 



var \ng = 4 [ 1 



1 i-i^ 

TT V A^ 



TT 



1- Vh+^(i) 



(15) 



For comparison, with diagonal disorder, (Ing) = —ADL and var Ing = SDL [Q. Although 
for even A^ the conductance distribution is log-normal both for diagonal and for off-diagonal 
disorder, there are some important differences: The presence of a term cx L^/^ in (In g) , the 
/9-dependence of (In (7), and the absence of universal fluctuations of Ing (they depend on 
rj) are special for off-diagonal disorder and do not appear in the standard case of diagonal 
disorder 0. 

What is rj for a particular microscopic model? Although we cannot answer this question 
in general, we can discuss some examples of microscopic models of quantum wires with hop- 
ping disorder and compare theoretical predictions from Eq. (|1^) with numerical simulations. 
We start from the Schrodinger equation for a set of coupled chains with random hopping, 
which in general can be written as 



n+1 



n— 1 • 



(16) 



Here n labels the position along the wire (in units of the lattice spacing), ipn is an A^- 
component wavevector, and t„ is an A^ x A^ hopping matrix, see Fig. ^ja. (This is a more 
general formulation than the one of Fig. |I]. Note that the lattice of Fig. |^ is in this class.) 
The hopping matrices tn are real (complex) for /3 = 1 (2). We connect the disordered region 
of length L to two ideal leads for n < and n > L, characterized by t„ = U. For zero energy 
it is possible to solve for the transfer matrix explicitly. 



M(L) 



m 



m 



L/2 
n=l 



t-l^ 



2n.-1^2n 



(17) 
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n=0 n=1 n=2 n=L n=0 n=1 n=2 n=3 

FIG. 2. (a) Most general random hopping chain corresponding to the Schrodinger Equation 
([T6| ) for = 3. (b) Random flux model: A square lattice where each plaquette is threaded by a 
random flux <I>. 



(For convenience, we assumed that L is even; recall that we use a basis of left and right 
movers.) The product m)m has eigenvalues exp(2a;j), j = 1, . . . , A^; the product m~^m~^^ 
has eigenvalues exp(— 2a;j). Being an arbitrary N x N matrix with real or complex elements, 
the matrix m is an element of the linear group GL(A^, M) or GL{N,C), and, taking into 
account the considerations outlined above, the L dependence of the transfer matrix M can 
be described by the random trajectory of the matrix m in GL{N, M) or GL{N, C). 

We now consider two examples. First, we consider the microscopic model used in Ref. 
|TI| (and also Ref. [^]). There, a distribution for the t„ was assumed that was invariant 
under unitary transformations of the chains. 



tr 



(18) 



Here Wn is a real (complex) matrix with independently and identically Gaussian distributed 
elements for /? = 1 (2). In this case, for small w, the distribution of the radial coordinates 
Xj of the transfer matrix M can be found explicitly. One obtains that the diffusion rates Dr 



are equal, and the parameters Xj 
and 7] 



w 



1 |0|,|4 



and Ds 
with Dch 

(|I^), including the dependence on 77, can be derived from Eq. 
a Gaussian distribution for the matrix Wn that involves correlations between the diagonal 
elements, 



obey a Fokker-Planck equation of the type (|T2D 
We mention that the most general Fokker-Planck equation 

using tn = expiWn) and 



(19) 



As a second example, we consider the random flux model [T^IT^ . This model corresponds 
to a square lattice, where each plaquette contains a random flux, see Fig. ^d. The hopping 
matrices t„ have the form 



/ 1 e'*^' 

1 
















1 e 



The phases 



bj^n are distributed such that the fluxes (^j^n 



l^j,n+l 



(20) 



bj^n are independently 



distributed. In order to find the distribution of the radial coordinates Xj for the random flux 
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W 



0.6 



0.8 



FIG. 3. The ratio c{N) = — linii^oo (In 5)/var In as a function of the disorder strength w of 
the longitudinal hopping parameters, for fixed randomness of the transversal hopping parameters. 
Statistical errors are smaller than the sizes of the symbols. 



model, we need to compute the matrix m, and find the eigenvalues exp(2a;j), j = 1, . . . ,N, 
of mm^, see Eq. (0). We do not need to solve this problem exactly: Each of the hopping 
matrices tn in the random fiux model has det t„ = 1, which implies exp(2A^a;) = det mm^ = 1 
for all lengths. Hence, the average radial coordinate x does not diffuse, so that = 0. We 
conclude that rj = for the random fiux model. 

The reason why = or rj = for the random fiux model is that det t„ = 1 for all n. 
More generally, all random hopping models which have det t„ = 1 are described by Eq. (|T^ 
with 7] = 0. One example of such a model is a random hopping model on a square lattice 
with only randomness in the transverse hopping amplitudes. In the general case, however, 
there will be both randomness in the transverse and in the longitudinal hopping amplitudes, 
so that T] > 0. 

To compare the Fokker-Planck equation (0) to numerical simulations, we considered 
the quantity 



c(iV) 



(In, 



f3N/2 



lim 

L~*oo var In g 



7V+(l-2/7r)(iV-2+2?7)' 



0, 



if N is even, 



if N is odd. 



(21) 



Here, we used Eqs. (|T4D and (^) for (In (7) and var \ng. In the standard case of on-site 
disorder, one has c{N) = 1/2 |0. We have compared Eq. ( pij ) to numerical simulations 
for a square lattice of width N and length L, attached to perfect leads. For (3=1, the 
transverse hopping amplitudes are taken from a uniform distribution in [—0.2, 0.2], whereas 
for f3 = 2, the transverse hopping amplitudes are complex numbers with amplitude uniformly 
distributed in [0, 0.1] and a random phase. In both cases the longitudinal hopping amplitudes 
are real numbers taken uniform in [1 — w,l + w]. We have numerically computed the ratio 
c{N) by taking an average over more than 2 x 10^ samples. Results for c{N) as a function 
of w for N = 2 and = 4 and for zero energy e = are shown in Fig. H. In the presence of 
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the chiral symmetry, i.e. for zero energy, we expect that c{N) depends on the details of the 
microscopic model (i.e., on the parameter w) through the parameter 77. For w = 0, we have 
?7 = 0, and hence c(2) = (3/2, c(4) = /?/(3 — 2/7r). As the randomness w in the longitudinal 
hopping amplitudes increases, we expect an increase of 77, and hence a decrease of c{N), 
which is confirmed by the numerical data shown in the figure. We have also shown data for 
energy e = 0.2 where the chiral symmetry is broken.^ In this case we find c{N) = 1/2 for 
all w, in agreement with the literature 0. [The slight increase of c{N) with w is attributed 
to a breakdown of the weak disorder condition.] 

Except for the case of the random flux model and the random-hopping model with 
transverse random hopping, where 17 = 0, we do not know how to compute the parameter 
7] explicitly. Numerical simulations show that, typically, the parameter 77 is of order unity. 
The effect of a nonzero t] on the conductance distribution is most pronounced for small A^. 
For large A^, the effect is small, because only the radial coordinates near zero contribute to 
transport and the additional correlations between these coordinates caused by a nonzero rj 
decrease as 1/A^. Therefore, in the limit of large A^, we still expect a universal conductance 



distribution. Here, we would like to remark that in Ref. ||2^, the Fokker-Planck equation 



(|T^) with ?7 = 1 was used to describe localization in the random-flux model. Despite the fact 
that the random flux model corresponds to 77 = 0, excellent agreement was found between 
the theory and numerical simulations performed for A^ = 15 and up. 

Before concluding, we would like to make two remarks. First, it is known that universality 
of one-parameter scaling breaks down in (quasi-) one dimensional disordered systems if there 
exist long-range correlations in the disorder, for instance in periodic-on- average systems ||2^ 



or for random-hopping chains in the presence of a staggering of the hopping parameter 
||10|| , as is relevant e.g. for narrow gap semiconductors or charge-density wave materials 
27| . (In the latter case, the staggering of the hopping parameter adds a drift term to 
the Brownian motion of x, whereas the Brownian motion of the relative positions of the 
Xj remains unaffected |T^.) The nonuniversality that we discuss in this paper is of quite 



a different origin: It follows directly from the geometric structure of the transfer matrix 
group; no long-range correlations in the disorder are involved, all the hopping amplitudes 
in the examples we considered have independent distributions. It should be noted that the 
appearance of a one-parameter family of scaling equations also occurs in the case of random 
Dirac fermions in two dimensions, where one finds a line of fixed points, rather than a single 



fixed point |llT|-|13 



Our second remark is of a more mathematical nature. Altland and Zirnbauer I12BIJ2B 



and Caselle |§] have argued that Cartan's classification of all symmetric spaces offers a 
complete classification of all possible random-matrix theories. These random-matrix theories 
appear in triplets: distributions of eigenvalues of Hermitian matrices (i.e., energy levels), 
distributions of eigenphases of unitary matrices (i.e., scattering phase shifts), and a Fokker- 
Planck equation for the radial eigenvalues of a transfer matrix. It is the latter kind of random- 



^ We have not shown the data points for w = and = 2 at finite energy, because in that 
case there is an extra reflection symmetry that is not taken into account in the standard DMPK 
equation. 
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matrix theories that we have considered here. While Cartan's classification is complete for 
all semi-simple transfer matrix groups, it does not take into account the phenomenon that 
we have presented in this paper: that for some disordered systems the transfer matrix group 
is not semi-simple, so that they cannot be represented by a single element from Cartan's 
table. 

To summarize, we have considered the Fokker-Planck equation for the transmission eigen- 
values of a quantum wire with off-diagonal (hopping) disorder from a geometric point of 
view. Under the same assumption of weak disorder that leads to the universal DMPK 
equation in the standard case of diagonal (on-site) disorder we have found that the 

Fokker-Planck equation for the transmission eigenvalues of a quantum wire with off-diagonal 
disorder contains an extra parameter that depends on the microscopic details of the disor- 
der. The existence of this extra parameter leads to a nonuniversality of transport properties 
in random-hopping chains, which is most prominent if the number N of coupled chains is 
small. 
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